Quantum anti-Zeno effect without rotating wave approximation 
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In this paper, we systematically study the spontaneous decay phenomenon of a two-level system 
under the influences of both its environment and continuous measurements. In order to clarify some 
well-established conclusions about the quantum Zeno effect (QZE) and the quantum anti-Zeno 
effect (QAZE), we do not use the rotating wave approximation (RWA) in obtaining an effective 
Hamiltonian. We examine various spectral distributions by making use of our present approach in 
comparison with other approaches. It is found that with respect to a bare excited state even without 
the RWA, the QAZE can still happen for some cases, e.g., the interacting spectra of hydrogen. But 
for a physical excited state, which is a renormalized dressed state of the atomic state, the QAZE 
disappears and only the QZE remains. These discoveries inevitably show a transition from the QZE 
to the QAZE as the measurement interval changes. 

PACS numbers: 03.65.Xp, 03.65.Yz, 42.50.Ct 



I. INTRODUCTION 

The quantum Zeno effect (QZE) is vividly described as 
a term "a watched pot never boils" in some quantum ver- 
sion [l[ . Usually it is used for describing a class of effects 
in which constant monitoring of a quantum system dras- 
tically slows down its dynamic evolution This may 
be a coherent transition (e.g., the Rabi oscillation Q) 
and an irreversible process as well. For instance, any un- 
stable state can be prevented from decay when adequate 
measurements are frequently applied to the system j^- 
7]. Here, the couplings to a reservoir would induce an 
exponential decay if there were no measurements. 

On the other hand, it was predicted that the decay 
could also be enhanced by frequent measurements ob- 
served under somewhat different conditions, leading to 
the so-called quantum anti-Zeno effect (QAZE) 
When the coupling to a surrounding environment (a 
reservoir) is taken into consideration, the generic QZE 
may not be attainable since the required measurement 
interval is out of reach in experiments. Furthermore, 
under the influence of the reservoir with some spectral 
distribution, the decay process could be significantly ac- 
celerated by continuous measurements. 

Recently it was recognized [l(| HH that the theoreti- 
cal prediction [8] for the reservoir-enhanced decay phe- 
nomena may be based on the rotating wave approxima- 
tion (RWA) 0, where the counter-rotating terms are ne- 
glected as they are high-frequency oscillating. A quite 
natural question follows as whether or not the existence 
of the QAZE really relies on the counter-rotating term, 
which is usually ignored in many applications since it 
possesses high frequency oscillation in the interaction pic- 
ture. 

In this paper, we will generally tackle this problem by 
investigating the influence of the counter-rotating term 
on the QAZE. Without making the RWA, as done in 
Ref. [10j, we develop a direct canonical transformation 



approach [12|, [l3[ to derive an effective Hamiltonian. It 
is equivalent to the second order perturbation approach. 
The obtained effective Hamiltonian is exactly solvable 
since it possesses the same form as that for the case 
with the RWA. Our calculation properly shows that for 
the spontaneous decay there exists a transition from the 
QZE to the QAZE as the measurement interval changes. 
In other words, with respect to the bare excited state 
(the product state of the atomic excited state and the 
vacuum of the reservoir) in the spontaneous decay, the 
counter-rotating terms are irrelevant to the occurrence of 
the QAZE for some spectral structures. As predicted, the 
essential difference between these approaches with and 
without the RWA could be disregarded in some cases. 

In addition to the spectra of hydrogen atom, we extend 
our research to the general situations with different kinds 
of spectral structures. Our finding shows that the QAZE 
seems to be universal except when some certain require- 
ment is met for a sub-Ohmic spectrum. Furthermore, in 
order to compare with the existing research [l(| , we also 
start from the same unitary transformation, but choose 
the bare excited state, which is different from the physi- 
cal excited state (the one excited from the ground state 
of the original Hamiltonian) in Ref. [10], as the initial 
state. Then the QAZE is again witnessed for the cases 
of hydrogen's spectral structure and others as well. The 
discrepancy between our result and the former one [lOj 
is attributed to the different choices of the initial states. 

The paper is structured as follows. In the next sec- 
tion, with a special transformation, we obtain the ef- 
fective Hamiltonian and thus the modification of the 
atomic spontaneous decay rate due to the counter- 
rotating terms. Sec. Ill discusses the transition from the 
QAZE to the QZE for different spectra. In Sec. IV, with 
the same initial state, we start from anther transforma- 
tion and arrive at the same conclusion for the hydrogen 
atom as the one in the previous section. And a brief sum- 
mary is concluded in Sec. V. Furthermore, we prove that 
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the survival probability of the atom in the excited state is 
equivalent to the survival probability of the initial state 
for the spontaneous decay in Appendix [A] In addition 
to Sec. IV, Appendix [B] presents the details about the 
calculation of the survival probability. 



II. EFFECTIVE HAMILTONIAN WITHOUT 
ROTATING WAVE APPROXIMATION 

We generally consider the QAZE for a two-level atom 
interacting with a reservoir in vacuum in the weak cou- 
pling limit. According to A. O. Caldeira and A. J. 
Leggett [lH, the reservoir weakly coupled to an open 
quantum system can universally be modeled as a collec- 
tion of many harmonic oscillators with annihilation (cre- 
ation) operator b k (b\) for fcth mode with frequency ui k . 
Let <J x ,y,z be the Pauli operators. And 



a = -(a x ±ia y ) 

are the raising and lowering operators for the atom with 
the excited state |e), the ground state |<?) and the energy 
level spacing Q, respectively. Then the total system is 
described by the Hamiltonian H = Hq + Hj: 



H Q = y^Ukbtbk + — O z , 



k 



Hi=J2sk[(b k + b{)a" 



h.c. 



(1) 
(2) 



Here, we have assumed the coupling constants g k s to 
be real for simplicity. However, we would like to say that 
the main result does not change if we start from a general 
assumption that g k s are complex numbers. 

As the interaction term Hi contains the counter- 
rotating terms, i.e., the high-frequency terms with fre- 
quencies ±(uj k + Q) like 

V = bl<7 + e* {u}k+n ^ + h.c. 

in the interaction picture, the Hamiltonian H is not ex- 
actly solvable even for the simple cases of single mode or 
single excitation. We use the generalized versio n |l2| of 
the Frohlich-Nakajima transformation exp(— S) [15|, [l6| 
to eliminate the high-frequency terms in the effective 
Hamiltonian. Here, 



(3) 



is an anti-Hermitian operator, where A^s remain to be 
determined. Up to the second order, the effective Hamil- 
tonian H c g = exp(— S)H exp(S') is given as 



H eS = H + H x + ^[H u S] + \[Hi,S] 



(4) 



Now we require b\,a + + h.c. to be eliminated from the 
first order term 

Ht = fli + [flo,S] 

][9k(b k + b\)a + + A k {uj k + n)b\a + + h.c] 



E 

*; 

E 

k 



gk(b k a + + h.c). 



The above equation gives the coefficients 

9k 



A, = 



U! k + il 



(5) 



Note that for a state \^>) which fulfills the Schrodinger 
equation before the transformation, i.e., 



H If) 



(6) 



we can prove that the state after the transformation 
I*) = exp(— S) |\&) satisfies the Schrodinger equation, 

(7) 



with the effective Hamiltonian 



H, 



cff : 



Hi 



u k b\b k + — cr z + 2_^gk{bkO + + h.c), (8) 

k 



where we have omitted the high-frequency intercrossing 
terms such as b\bu, and b k /b k , and the modified level 
spacing for the atom is 

(9) 



ni = n + Y-°±-. 
Y uj k + n 



Here, the shift of the level spacing can be regarded as the 
Lamb shift, and also called the AC stark modification in 
atomic physics and quantum optics. Furthermore, in the 
above calculation, the term 



is replaced by 



E 



9l 



LJ k + 



since for the single excitation case its contribution results 
in a small modification in the Vih mode gf/(oJi + f2). 
We remark that, for those modes k ^ k' with smaller 
frequency differences, the terms b\b k i (k ^ k') could have 
larger contributions in quantum dynamics, but for some 
initial states we will choose them to be of the second 
order. This problem has been considered in Ref. (loj . 

Let us first consider the QAZE for the spontaneous 
decay process where the initial state can be chosen as 
l e ; {0}) — l e ) ® |{0}) w hh the atom in the excited 
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state \e) and all modes of fields in the vacuum state 
|{0}) = Y\k ® -Due to the special unitary transfor- 
mation exp(— S), the initial states before and after the 
transformation are identical, i.e., 



{0}) = (I 



S+^ 2 )|e,{0}) = |e,{0}). (10) 

We note that, for the generalized Frohlich-Nakajima 
transformation in Ref. jToj . the initial state would be 
changed. For other cases, it will be illustrated that the 
uses of changed and unchanged initial states would result 
in the different conclusions about the discussions of the 
QAZE. 

When the atom is projected onto the excited state pro- 
vided that the total system evolves from the initial state 
e, {0}), the survival probability is 



P(t) = \x(t)\' = Tr(\e) (e\e-^\ 



,{0})(e,{0}\e mt ). 

Thus, as shown in Appendix |Al the survival probability 
after n measurements 

( e ,{0}| e -^|e,{0})| 2 " 

{e,{0}\e s e-^^-^ 



P(t 



(e,{0}|e 

-Rt 



-Wear 



■{0}) 



e,{0}}| 

2n 



(12) 

is calculated in the transformed representation where the 
new initial state just coincides with the original one. 
Here, the decay rate [|| 



R = 2tt 



du}F(u},fti)G(u) 



(13) 



is the overlap integral of the measurement-induced 
atomic level broadening 



1 f 
F(w, f2i) = — / 
2tt J_ t 



-smc 



dt(i-^y 2nit 9(T-\t\y 

T 

(lu — ^i)r 



2tt ' L 2 
and the interacting spectrum 



(14) 



dt y £k e i(u-u k )t = 



E- 2 

A,- 



g k S(u}-u} k ). (15) 



The above obtained result seems to be the same as that 
in Ref. Q, but the essential difference is that the peak 
of F(lj, Oi) has been shifted due to the counter-rotating 
terms. In this approach for practical problems, this shift 
may not have significant effect on the physical result (see 
the illustration in Fig. In the physical systems that 
we are considering, i.e., hydrogen atom, the influence of 
the counter-rotating terms is tiny small since the energy 
shift |Oi — fi| is relative small with respect to the dis- 
tance between the original energy level spacing Q and 
the peak of the interacting spectrum ujq- However, there 
may appear some different results in artificial systems 
such as circuit QED [l?} • We will check this observation 
for various cases as follows. 




FIG. 1: (color online) Schematic of the difference between 
the overlap integrals with and without the RWA. The blue 
solid line for the measurement function F(u), O) centered at 
the original atomic level spacing Q. for the case with the RWA 
and the red dotted line for the measurement function F(uj, £7i ) 
centered at the modified frequency fli for our current case 
without the RWA, the green dashed line for interacting spec- 
trum G(u>) centered at u>q. 



III. QUANTUM ANTI-ZENO EFFECT FOR 
DIFFERENT INTERACTING SPECTRA 

Having obtained the effective decay rate modified by 
the counter-rotating terms, we examine the above obser- 
vation for specific spectra in investigating the QZE and 
the QAZE. 



A. Quantum Anti-Zeno Effect for Hydrogen Atom 

Let us first investigate the decay rate for the hydro- 
gen atom in the vacuum of electromagnetic fields. We 
consider two usual transitions of the hydrogen atom, i.e., 
2p-ls and 3p-ls, with the interacting spectra (Til. [l9j 



G 



2p-ls 



[l + (^) 2 ] 4 ' 



and 



rA4l + 2(^)T 

respectively, where 

r\ = 6.435 * 1CT 9 , uj c = 8.491 * 10 18 rad/s 
r; ' = 1.455 * 1CT 9 , J c = 7.547 * 10 18 rad/s 



(16) 



(17) 



(18) 
(19) 



The numerical calculations of the decay rate in Eq. (fT3")) 
are shown in Figs. [2] and [3l Here, we observe the emer- 
gence of the QZE and the QAZE as well. Starting from 
a large enough value of r, as the measurement interval 
decreases, the decay rate will experience an ascending 
procedure at the first stage. Since the decay rate is big- 
ger than the unperturbed one 

R = 27^(00, 



4 




since F(uj,SY) has already covered the main part of G(u>). 
And afterwards the opposite phenomenon is observed. In 
these two figures, also shown arc the differences between 
the decay rates obtained from the one with the RWA 
and the one without the RWA, AR = \R— i? rwa |. Notice 
that AR's are of the order of 10 -6 (in units of Rq). It's 
a reasonable result since the only effect of the counter- 
rotating terms lies in the modified level spacing And 
the small correction is of the order of 10 ~ 6 with respect 
to the original level spacing. 



FIG. 2: (color online) The decay rate vs measurement in- 
terval r for the 2p-ls transition of the hydrogen atom. Here, 
blue solid line for R and red dashed line for AT? = \R — R Ivla ,\. 
[fii - = 1.71 x 1(T 6 , fi = 1.55 x 10 1B rad/s and 

uj c /Q = 550. Notice that AR is enlarged by 10 6 times. Notice 
that AR is enlarged by 10 6 times. In all figures, the measure- 
ment interval r is in units of atomic level spacing I/O. and the 
decay rate R is normalized with respect to the unperturbed 
one Ro. 



B. Quantum (Anti-)Zeno Effect for General 
Spectral Distribution 

Afterwards, we generally investigate the QAZE for dif- 
ferent spectral structures. Especially, we discover the 
condition when the QAZE disappears. In general, the in- 
teracting spectra are classified as three categories. They 
can be written with a uniform spectrum function [20| 



(20) 



where A is a constant and w c is the cutoff frequency. 
For an Ohmic spectrum, s = 1 while s < 1 and s > 1 
for sub-Ohmic and super Ohmic spectra respectively. In 
Fig- HI the transition from the QAZE to the QZE is again 
observed. It is a predictable result since the peak of the 
spectrum function is located at w = su) c . As long as 
soj c 3> £1, the QAZE definitely occurs. Moreover, on 
condition that 



0.001 I 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

0.0 0.2 0.4 0.6 0.S 1.0 



FIG. 3: (color online) The decay rate vs measurement in- 
terval r for the 3p-ls transition of the hydrogen atom. Here, 
blue solid line for R and red dashed line for AR = \R — R Ivla \. 
|Oi - n\/n = 1.27 x 1(T 6 , n = 1.83 x 10 16 rad/s and 
ui' c /Q = 412. Notice that AR is enlarged by 10 6 times. 



the QAZE occurs before it reaches the climax. After the 
turning point, the trend is changed. It is obvious that 
the decay rate drops steeply as the r is further reduced. 
When the normalized decay rate falls below 1, the QZE 
is present. As the measurement becomes more and more 
frequent, i.e., r — > 0, we observe the transition from the 
QAZE to the QZE. Mathematically speaking, the decay 
rate is the overlap integral of the measurement-induced 
atomic level broadening F(u>, fi) and the interacting spec- 
trum G(uj). F(lu, Q) is peaked at Q with width 1/r while 
G(lo) is maximized at a frequency of the order of the 
cutoff frequency uj c which is much bigger than the atomic 
level spacing ft. As r decreases from a large enough value, 
F(uj, f2) covers more and more raising part of G(lo). As a 
consequence, the decay rate is enhanced and the QAZE is 
witnessed. When the measurement interval r is reduced 
to the order of 1 /u) c , the decay rate will no longer increase 



fii ~ tt = suj c , (21) 

the QAZE is wiped out and only the QZE takes place, 
as shown by the black dot-dashed line in Fig. Addi- 
tionally, the difference between the decay rates with and 
without the RWA is plotted in Fig. [5] Since the contri- 
bution is no more than 10~ 3 for the given parameters 
A = 10~ 8 and oj c /CI = 500, the counter-rotating terms 
thus can be neglected as the routine work done in quan- 
tum optics. 

However, the above analysis is based on the assump- 
tion of a small energy level shift. For some physical sys- 
tems, this shift may play an important role in the ex- 
istence of the QAZE. For a given interacting spectrum 
as shown in Eq. (|20p . the modified energy level spacing 
reads 

fij = Q + Ae£r2 s r(l + s,0)r(-s,— ), (22) 

where 

/>oo 

T(u,z) = / t^e^dt. (23) 

J z 

As stated above, the QAZE disappears if the peaks of the 
measurement-induced atomic level broadening and the 
interacting spectrum coincide. Therefore, the distance 



FIG. 4: (color online) The decay rate vs measurement in- 
terval r for different spectra. Here, (a) green dotted line for 
super-Ohmic spectrum s = 2, (b) blue solid line for Ohmic 
spectrum s = 1, (c) red dashed line for sub-Ohmic spectrum 
s = 0.5 and (d) black dot-dashed line for sub-Ohmic spectrum 
s = Q,/u c . For all spectra we set A = 10 -8 , U) c /Cl — 500. 



FIG. 6: (color online) The frequency distance AQ — Qi —slu c 
vs the parameters A and s with u c /£l = 500. 



IV. ALTERNATIVE APPROACH FOR THE 
DECAY PHENOMENON OF HYDROGEN ATOM 




In this section, with a different choice of the transfor- 
mation exp(— S'), we will obtain the same result as the 
preceding section for the hydrogen atom. This transfor- 
mation approach can also work well for other spectral 
distributions but here we do not repeat the straightfor- 
ward calculations. 

We choose the same transformation U' = exp(— S') as 
that in Ref. [Io[ with 

S' = Y,M(bU + ~ b k a-) + (bl<7- - b k a + )], (24) 



FIG. 5: (color online) The decay rate difference AR = \R — 
i?rwa| vs measurement interval r for different spectra. Here, 
(a) black dot-dashed line for sub-Ohmic spectrum s = Q/u> c 
with |fii - f2|/fi = 1.55 x 10~ 3 , (b) red dashed line for sub- 
Ohmic spectrum s = 0.5 with |fii — = 9.41 x 10~ 4 , 
(c) blue solid line for Ohmic spectrum s — 1 with — 
= 6.57xl0" 4 , and (d) green dotted line for super-Ohmic 
spectrum s = 2 with \Qi — = 4.00 x 10~ 4 . Parameters 
are the same as those given in Fig. [4] 



AO = 51i — suj c between these two peaks is plotted vs 
the parameters A and s in Fig. [SJ As shown, the distance 
AJ7 increases monotonically with increasing A. This is 
because the energy level shift, i.e., the second term on 
the right hand side of Eq. (|2"2")l . is proportional to A. 
Physically speaking, the larger the parameter A is, the 
stronger the coupling between the atom and the reservoir 
becomes. As a result of the stronger coupling, the energy 
level shift is enlarged. Besides, it is seen that Aft falls 
as the parameter s raises. Thus, for a matching pair of 
A and s, the two peaks of -F(w,J7i) and the interacting 
spectrum G(w) are the same. In this case, there will be 
only the QZE. Besides, we also notice that the QZE was 
well explored for a two-level system in either a low- or 
high-frequency bath beyond the RWA [2l| . 



and 



'9k 



(25) 



to eliminate the counter-rotating terms b\.<T + + b k cr~ in 
the desired effective Hamiltonian. This transformation is 
different from ours in that it includes the slow-oscillating 
terms, i.e., b\a~ — b k a + . As a consequence, it modifies 
not only the atomic energy level spacing, but also its cou- 
pling to the reservoir and thus the interacting spectrum. 

By virtue of omitting higher order terms, i.e., b\.b\,, 
bkbk' , the effective Hamiltonian is straightforward given 
as 



H'a 



E 

it 



Lt.. ^ 
WfeO^Ofc + —<7 Z 



k 



h.c. 



where the modified coupling constant is 

2fl 



9k = 



:9k, 



LOk + O ' 

and the modified atomic energy level spacing is 

QgkA k 



ft' = + 2 ^ 



(26) 



(27) 



(28) 
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with the coefficients 



with the "renormalized" free Hamiltonian 



A k 



-9k 



We would like to mention that the term 



(29) 



EVlg k A k 



is replaced by 



E 



UJ k + CI 



in the above calculation since their contribution results 
in small modification. 

Since the original Hamiltonian H and a state \tp(t)) 
fulfill the Schrodinger equation H\i[;(t)) = idt\ifj(t)), we 
emphasize that it is the transformed state \ip(t)) s — 
exp(—S')\i/j(t)) and the effective Hamiltonian H' eS that 
meet the same requirement H' eS \if)(t)) s — idt\ip(t)) s . 
Thus, in general cases the initial state \ip(0)) before 
the transformation should be changed as |?/>(0)) s = 
exp(— S')\ip(0)) after the transformation. 

In practice, the choice of the initial state relies on the 
concrete physical problem. So far as the vacuum-induced 
spontaneous decay is concerned, we should choose the 
bare excited state. We emphasize that this choice is con- 
sistent with the one in Ref. We also remark that 
it would be more reasonable to start from the same ini- 
tial state when we refer to the influence of the counter- 
rotating terms on the QAZE. Besides, we can also choose 
the physical excited state elsewhere, i.e., in Ref. ficj . 
It is a reasonable consideration since the ground state 
\g, {0}) of the Hamiltonian under RWA is replaced 
by exp(S")|<7, {0}} due to the presence of the counter- 
rotating terms in the interaction Hamiltonian ([2]) (22j. 
Therefore, the initial state may be exp(S")|e, {0}) instead 
of e, {0}) under the condition that the initial state is 
prepared from the ground state exp(S')\g, {0}} through 
excitation by laser. These two different choices will result 
in distinct consequences. 

The problem is solved in the interaction picture with 
respect to 



Un =. 



iHit 



(30) 



K = ^u k b\b k + —a z . 



(31) 



And the interaction Hamiltonian H\ = UqH'-JJq = 
U (H' oS - H' )U+ reads 

H[ = ^2g' k {bk<? + e m '~ UJk)t + bla-e-W-"^*). (32) 

k 

The time evolution of the wavefunction 

W t {t)) = a(t) |e, {0}) + J2 &(*) |s, 1*> (33) 

k 

is governed by the effective Hamiltonian H[ , 

idt Wj) = Hi Wj) . (34) 

Here, \g, l k ) denotes the atom in the ground state \g) and 
one excitation in the fcth mode. 

Then, the coefficients meet the following demands 



p k = -ig' k ae- i{ - n '- Uk ^. 



(35) 
(36) 



However, when calculating the survival probability, we 
should return to the Schrodinger picture, i.e., 



— i — t 

a -> ae 2 , 



0k ^e-*( w *-T)*. 

We would like to remark that with the initial state 
e,{0}), the considered survival probability for the ex- 
cited state of the atom under the original Hamiltonian 
H is 



P(t) = \x(t)\< = Tr(|e) <e| e~* m \e, {0}) (e, {0}| e^) 



(37) 



Correspondingly, the effective Hamiltonian and the ini- 
tial state after the above transformation are H' eS and 
e~ s |e,{0}), respectively. Then, one has the survival 
probability amplitude 



x (t) = (e, {0}| e^V^e" 3 ' |e, {0}} 

~d (e, {0}| e- iH '^ \e, {0}> - C 2 £ A k (e, {0}| e^ H '^ \g, l k ) - C 2 £ A k (g, l k \ e^* 4 |e, {0}) 

k k 

+ J2Al{9,h\e- iH '°« t \g,lk), (38) 



where we have dropped the off-diagonal terms for the fourth term on the right hand side and 

Cx = (l-i]T4) 2 ^l-£4, (39) 



2 



C 2 = l-W A l ( 4 °) 
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In Appendix El we present the detailed derivation in ob- 
taining the first line of Eq. ([38]) from Eq. ([57]). After 
a series of deductions, the survival probability after n 
measurements is written as 



P{t = nr) = \x{r)\ 2n = e- R \ 
where the decay rate 

/>oo 

R = 2tt dujF(uj,n')G'(uj) 



(41) 



(42) 



is the overlap integral of the measurement-induced 
atomic level broadening 



J . 2 / w — ^' 



F(u,Q') = — sinc 2 ( 
and the interacting spectrum 

r;'u-i = X;/(«ik)5**(«- 



r) 



with 



/(w fc ) = 1 + 



(3ft-0' + 2w fc )(ft-ft') 

(wfe + fi) 2 



(43) 



(44) 



(45) 



Notice that the measurement-induced atomic level broad- 
ening differs from ours in Eq. (|14|) in that it is centered 
at a different modified level spacing 57'. And the inter- 
acting spectrum is also altered with an additional factor 
f(tu) in contrary to the unaltered one in Eq. (|15p . For 
the necessary details, please refer to Appendix iBl 

For the hydrogen atom, the modified interacting spec- 
trum of the 2p-ls transition is 



r/ui 



214 



(46) 



with 77 and uj c already given in Eq. (ITS)) . Notice that the 
second term in Eq. (|45)) is a small correction to the one 
with the RWA of the order of 10 -8 . This is in consistence 
with the orders of the numerical results of AR in Fig. 05J 
The relation between the decay rate R and the mea- 
surement interval r is plotted in Fig. [7) As shown in this 
figure, the decay rate is well separated into two parts 
with the climax being the boundary. In the left part, as 
the measurements are done more and more frequently, 
i.e., r — > 0, the decay rate falls monotonously. When it 
is less than the decay rate without measurement 

R' Q = 27rG / (0 / ), 

i.e., R/R'q < 1, the QZE takes place. To the right of the 
climax, the opposite trend is witnessed. In this region, 
the shorter the measurement interval is, the larger the 
decay rate is. Since R/R'q > 1 for the whole region to the 
right of the climax, one can observe the QAZE, which was 
predicted to be obliterated due to the counter-rotating 
terms 10]. The discrepancy between their result and 




FIG. 7: (color online) The decay rate vs measurement in- 
terval t for the 2p-ls transition of the hydrogen atom. Here, 
blue solid line for R and red dashed line for AT? = \R — T? r wa|. 
\Q! - Q\/Q = 5.69 x 1(T 8 and lu c /Q, = 550. Notice that AR 
is enlarged by 10 s times. 




FIG. 8: (color online) The decay rate vs measurement in- 
terval r for the 3p-ls transition of the hydrogen atom. Here, 
blue solid line for R, red dashed line for AR = \R — R Iwa .\ and 
green dotted line for the result from Ref. 
5.32 x 10" 8 and uj c /Q, = 412. Notice that AT? is enlarged by 
10 s times. 



ours is attributed to the different choices of the initial 
states [10| . Here, we also emphasize that our approach 
is very simple and concise in contrast to theirs, which is 
shown in Appendix [Bj 

For more evidence, we resort to the 3p-ls transition of 
the hydrogen atom. The interacting spectrum is adjusted 
as 



212 



G' 3p _ ls (u) = /(«) 



[i + (-5) 



210 



(47) 



where r{ and uj' c are given in Eq. (|19|) . The result of 
the numerical calculations is displayed in Fig. [5] Here, 
we again observe the complete opposite predictions of 
the QAZE in contrary to Ref. [13]. Moreover, the RWA 
offers a good approximation in the weak-coupling limit 
since the disagreement between its and the exact result 
is trivial. 



8 



V. CONCLUSION 

In summary, we have studied the role of the counter- 
rotating terms of the atomic couplings to the reservoir in 
the irreversible atomic transition during the continuous 
measurements. By the generalized Frohlich-Nakajima 
transformation, the exactly solvable Hamiltonian is ef- 
fectively obtained without the RWA in the form of the 
standard "spin-boson" model. We discovered that when 
we consider the spontaneous decay of the bare excited 
state even without RWA, the QAZE remains if the proper 
measurement interval is given. And a transition from the 
QZE to the QAZE as the measurement interval changes. 
As for the findings in Ref. [10( and its following papers, 
it is observed that the disappearance of the QAZE under 
the approach without RWA is mainly due to the choice of 
the physical excited state. In Ref. [l(| , this initial state is 
the excited state of the renormalized Hamiltonian, which 
is essentially an entangled state of photons and atomic 
states. For the physical systems in realistic world, the 
influence of the weak-coupling counter-rotating terms on 
the decay rate is tiny small and can be negligible. We 
have utilized two different approaches for the generaliza- 
tion of the Frohlich-Nakajima transformation. For the 
same initial state, i.e., the bare excited state, consistent 
conclusion is obtained for the interacting spectra of the 
hydrogen atom. By comparing the two effective Hamil- 
tonians, we find out that in our approach there are only 
one parameter modified in contrary to one more set in 
Ref. [ioj . To further verify the universality of the pres- 
ence of the QAZE, we also extend our investigation to 
different types of spectra. It is discovered that when the 
cutoff frequency and the atomic level spacing fulfill some 
condition, only will the QZE emerge for the sub-Ohmic 
spectrum. We notice that by means of the QZE in the 
dynamic version a quantum switch was proposed to con- 
trol the transport of a single photon in a one-dimensional 
waveguide under the RWA [231. However, we might look 
forward to some novel features if no RWA is invoked. 

Besides, it is worth underlining that the choice of the 
different initial states depends on the specific physical 
problem. So far as the QAZE for the vacuum-induced 
spontaneous decay is concerned, we should choose the 
bare excited state in that it is the vacuum that induces 
the spontaneous decay of the atomic excitation. On the 
other hand, on account of the preparation of the initial 
state, the physical excited state might be a better choice 
as well because it can be feasibly excited from the ground 
state of the original Hamiltonian. 

This work was supported by NSFC through grants 
10974209 and 10935010 and by the National 973 program 
(Grant No. 2006CB921205). 



Appendix A: Survival Probability P(t) 

When we refer to the QZE and QAZE, we take a 
projective measurement on the atom. Thus, we shall 
trace over all the possible states of the fields. Based on 
the above considerations, we give the detailed deduction 
about the survival probability P(t) — \x(t)\ 2 . 

Before the transformation, the original Hamiltonian is 

H = J2 u k b\b k + ^a z + J2 9k(b k + b\){a+ + a~) (Al) 
fe fc 

with the chosen initial state to be 

|*(0)) = |e,{0}). (A2) 

Then, we take a unitary transformation e~ s with 

S = J2 A k(bU + ~b k a-), (A3) 
fe 

the Hamiltonian is approximated to the second order as 

H ca = e- s He s 

= uj k b\b k + —a z + g k (b k a + + h.c), 

k fe 

(A4) 

in company with a transformed initial state 

|*(0)) S = e- s |*(0)} 
= e- s |e,{0}) 

= |e,{0}>. (A5) 

As a result, the evolution of the state reads 

|*(i)) S = e- iH ' tft |*(0)) S . (A6) 

When calculating the survival probability for the atom 
in the excited state, we shall return to the original picture 
and the density matrix for the total system is straight- 
forward given as 

P {t) = m))m)\ 

= e s \V(t)) SS (V(t)\e- S 

= e s e-' lH " !it e- s |*(0)) (*(0)| eV^e^ 

= e s e- iH " tft e- s |e, {0}) (e, {0}| e s e lH ^ l e~ s . 

(A7) 

The reduced density matrix for the atom is traced over 
the degrees of fields, i.e., 
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p s (t) = Tr B |*(t)) (*(£)| 

= ({0}| e s e-^ ff *e- s |e, {0}) (e, {0}| e s e iHMt e~ s |{0}) 
+ ^2(lk\ e s e- lH ^ l e- s |e, {0}) (e, {0}| e 5 e l ^«'*e" s |l fc ) 



S -iH B «t e -S | 6j { }^ |o}| e S e iH e „t e -S | lfclfc/ ^ 
I 



Therefore, the survival probability of the excited state of the atom is 
I 



(A8) 



P T(t) = Tr s (\e)(e\p s (t)) 



<e, {0}| e s e- m ^e- s |e, {0}) <e, {0}| e s e m ^e~ s |e, {0}} 



^S-iH aii t -S 



3,{0})(e,{0}|e 



+ £ (e, 1*1* I e s e- l ^ ft e- s |e, {0}) (e, {0}| e s e lH ^e~ s \e, l fc l*> 

I 



(A9) 



In the following deductions, we will show that multi- Then, the third term on the right hand side of Eq. (IA9 
pic excitation terms can be omitted as they will lead to equals 



small corrections to the final result. For the case of two 
excitations, 

e- s \e,l k ) ~ (I-S+^S 2 )\e,lk) 

= |e, l fe ) + A k \g, {0}) AkAk ' l e ' lfe ') ' 

v 

(A10) 

is a superposition of states with the total excitation of 
an even number, while 



e s \e, 



{0}> ~ (/-5+|s 2 )| e ,{0}> 

= |e,{0}) 



(All) 



has only one excitation. On account of -ffeff's property 
of conserving the total number of excitation, the second 
term in Eq. (|A9I) vanishes. 

For the case of three excitations, 

e~ s |e,lfelfc'> 
* (I-S+^S 2 )\e,l k l k ,) 
= |e, lfelfe/) + A k \g, l k ,} + A k , \g, l k ) 

~\ Y^( AkAk " l e ' + A k>A k " \e, lfelfe"))- 

k" 

(A12) 



J2\(e,lklk'\e s e- lH ^e- s \e,{0}}\ 



k.k> 



k'\e 



-iH a{t t 



e,{0})| 



k,k' 



+2 



^2A k {g,l k \e- iH ^\e,{0}) 



(A13) 



In the interaction picture, the interaction Hamiltonian 
reads 

H 1 = ^9k(bk<? + e l{ni ^ k)t + b\a-e- l{n ^ k)t ). (A14) 



The time evolution of the wavefunction 

|^> = a(t) \e, {0}) - £ AW |3, U> • (A15) 

fc 

is governed by the Hamiltonian H 1 , 

i9 t |^> = ff' ■ (A16) 

Straightforward, we attain the equations for the coef- 
ficients as 



i(ni-a;fc)t 



(A17) 
(A18) 
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The first term of Eq. (|A13|) is equivalent to 
%J2k k' ^1 \h'{t)\ when a(0) = 1. We can formally in- 
tegrate Eq. (|A18|) and replace a(t') by 1 to have 



o 



l l dt'g k e 
lo 

e -i(fl 1 -uj k )t _ ^ 



i(fii-u fc )t' 



9k - 



fii - LO k 



-2 sin 



2 (Qi-u k )t 



9k- 



— i sin(Oi — Lo k )t 



ill - Wfe 

As a result, the first term of Eq. (|A13|) 
2Y,Al\p k ,{t)\ 2 



(A19) 



k,k' 



= *E4E^ 



4 sin 



4 (ni-^,)t 



+ sin 2 (f2i — ujb')t 



(Qi - LO k >) 2 



pao) 



can be neglected based on the following considerations. 
For one thing, the first term in the second summation is 
proportional to i 4 in the short time limit and thus can be 
omitted. For another, the factor J2 k A\ is a small quan- 
tity, i.e., typically of the order of 10 -8 for the hydrogen 
atom. 

Moreover, the second term on the right hand side of 
Eq. (IA13[) is equivalent to 2 A k j3 k {t) exp(— iu> k t)\ 2 
where a(0) = 1 and the factor exp(— iLo k t) is due to trans- 
formation back to the Schrodinger picture. Here, 



k 

2 sin 2 ( " 1 ~" fc)t - i sin(Oi - w k )t 



k 

4 />oo 



= — [ / dusmc 
2 



Ql - LO k 
2 (Oi — Lo)t s—^ VL\ — LO 2 



E 



fil + LO 



g k S(Lo~Lo k )Y 



where the first term is proportional to t 4 in the short 
time limit, and the second term is of higher order with 
respect to the first term in Eq. (IA9|) . 



As a consequence, the second and third terms on the 
right hand side of Eq. (|A9|) can be neglected. In other 
words, the contributions from the multiple-excitation 
terms result in a small correction to the final result and 
thus we have 



p?(t) = \(e,{0}\ 



■{0}>| 



(A21) 



which is exactly the same as the one in Eq. (TT2"j) for a 
single measurement. 

For the second approach in Sec. IV with 



S' = Y j A k [{b\cj+ - b k a~) + {b\a- - b k a+)}, (A22) 



k 



+2t 2 \ [ diosmct^ - u)iV — ^] 2 we still have Eq. K§ 



P T(t) = Tr s (\e)(e\p s (t)) 

~ (e, {0}| e s 'e-^*e- s ' |e, {0}) (e, {0}| e s V^«' e - s ' |e, {0}) 
+ E <e, U| eS'e-^e- 8 ' \e, {0}) (e, {0}| e^e^e" 5 ' \e, l k ) 

k 

+ E (e, 1*1* I e s V^«* e - s ' ^ {Q}) ^ {Q} | e s' e ^ e - s ' ^ lfelfc;) + . . . _ (A23) 

k,k> 

I 

In this case, the effective Hamiltonian is replaced by while the transformed initial state is given as 

H' eS = e- s 'He s ' |*(0)} S ' = e~ S ' |tf(0)> 

= J2^ b l b » + Y az+ T,9'k(bkV + + h.c.), = e " S > ; {0}) 

* [/-5' + ^(^) 2 ]|e,{0}) 
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Then, wc obtain the transformed state at time t, 



As 



e" S '|e,l fc ) 



\^>'{t)) S ' = e- lH '^ |*'(0)} 5 ' 



[/-5' + -(5') 2 ]|e,l fc ) 



(A26) 



fc' fc 

-^Ay \g, lfclfe/) -^AfeAfe' |e,l fc '> 
fc fc' 

+ | ^ Afe/Afc. |e, lfclfe'U") (A27) 



and 



e- s '|e,{0}) ~ [I-5' + i(5') 2 ]|e,{0}) 



= (l-i^A 2 ,)|e,0)-^A fc | ff ,l fc ) 
fc' fc 

+ i^A fc ,A fc |e,l fc lfc') (A28) 



are of even and odd numbers of total excitations respec- 
tively, we have a vanishing second term on the right hand 
side of Eq. (1A"23)) . 

Because the contribution from 

e~ s ' e, lfclfc') 
~ {I-S' + ±(S') 2 ]\e,lklk>) 
= |e, l fc lfc') + A k \g, l fc ,) + A w \g, l k ) 

-^A k „ |5,lfclfc'lfc»> - ]^^A k A k ,, |e, lfc/ lfc//) 
fc" fc'/ 

|e,0) + i ^ Afc//vlfc/// |e,lfclfc/lfc//l fc ///) 
fc",fe/// 

-)^A k A k „ |e,lfe/l fe //) - ^J2 A k' A k" 
fc" fc" 

-i^A fc /A fc » |e,l fe lfc") - |e,lfclfc') 
fc" fc/' 

(A29) 



is of the order higher than A 2 ,, we can further omit the 
third term of Eq. f&23]) . 

In total, we can neglect all nonzero-photon terms in 
Eq. (|A23[) . Thus, based on the above calculations, we 
still have 



p?(t) = \x(t)\ 2 = (e,{0}| 



3,{0}> 



(A30) 



for the second approach, where the survival probability 
x(t) is the same as Eq. (|3"5]) . 

Judging from the above reductions, we may safely ar- 
rive at the conclusion that for both cases the survival 
probability of the atom in the excited state coincides with 
the survival probability of the initial state |e, {0}) and the 
single excitation approximation is reasonable. 



Appendix B: Survival Probability Amplitude x(t) for 
the Second Approach 



As shown in Eq. (|38|) . the survival probability x(t) 
is the summation of four terms due the modified initial 
state after the transformation. Since the calculation of 
the first term in Eq. (|38[) was already shown elsewhere, 
i.e., Ref. Q, we offer the detailed calculation of the re- 
maining parts in addition to the first term. 

For convenience, we multiply x(t) [Eq. by a factor 
exp(if2't/2), and have 



,(t)e 1 - 1 =C X (e, {0}| e- m ^ |e, {0}) e l ~ l - C 2 ]T A k (e, {0}| e"*^' \g, lfc) e^ 4 

fc 

- C 2 Y,A k (9, lfc I e"^ fft |e, {0}) e^* + £ A\ (g, l fc | e~ iH ^ \g, lfc) e^ 4 . (Bl) 



On the right hand side of the above equation, the first term 



(e, Wle-^'le, {0}) 



(B2) 
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is equivalent to a(t) when a(0) = 1. Eq. (j36|) can be 
formally integrated to yield 



Jo 

By substituting it into Eq. (|35[) . we have 



*=-E fti'ttk) 



2 ae -i(n'-« fc )t' e »(n'-ufc)t < 



(B3) 



(B4) 



For a sufficient short time i, we can replace a (£2) with 
a(0) = 1 and thus 



Jo Jo k 
= i-4(t). 



(B5) 



When transforming it back to the Schrbdinger picture, 
we have 



(e,{0}|e^« t |e,{0})=a e - 



(B6) 



By multiplying a factor exp(iri'i/2), the time-dependent 
factor is canceled, 



(e, {0}| e- lH ^ \e, {0}> e l ~ l = a. (B7) 

And the following function will be used in the calculation 
of x(t), 



^2 -iui k t' ' j 



t 9 ,u> — f2' 
2^ SmC ( ^ ' } < 



2ReJ„(t) = 2Re[t / dt'(l - ^)e* n '<' 

/oo 
dwF(u;,fi') G i( w )> (B8) 
-00 

which is the overlap integral of the measurement function 

1 r°° \t'\ , , , 

f{lj, n') = — dt'(i - u. * e(t - \t'\) e - lult 
27r J -00 1 

(B9) 

and the interacting spectrum 

G 1 (u) = — dt'^kfe-*"** e*** 
Z7r J-00 k 

= J2(9'k?5("-"k)- (BIO) 

He 

Before calculating the second and third terms in 
Eq. (|B1I) , we prove these two terms to be equal to simplify 
the calculations. For a general Hamiltonian H, which is 
time-independent and satisfies H ab = (a \ H \ b) = H ba for 



any two states \a) and \b) in the complete Hilbert space, 
one has 



F ab = (a | e~ im \b) 

{-itH) 7 



WE" 



\b) 



t ^(-itH) 2n |n . |V , (-itH) 2n+1 „ 



(2n)! 

= Re(F o6 ) + Hm(F ab ). 



(2n + l)\ 



(Bll) 



It is obvious that the first term on the right hand side of 
the above equation is real and the second term is pure 
image. On the other hand, one has 



,— iHt 



(b\e — |o, 
\(b\e- iHt \a)y 
l(a\e im \b)? 



n 

HE 



(itH) 



2n 



(itH) 



2n+l 



(2n)\ " (2» + l)l 

(-itH) 2n ^ (-itHf 



\b) 



(2n)\ 

= [Re(F ab ) - ilm(F ab )}^ 
=Re(F ab ) + i!m(F ab ) 
=F ab . 



(2n+ 1)! 



\b) 



(B12) 



The above condition H ab = H ba (for any two states \a) 
and I b)) for a general Hamiltonian H, means g% — gT in 
our current case. In short, one has 



(e, {0}| e-^«* | 5 , l k ) = (.g, l fe | e"^ °«* |e, {0}) . (B13) 



For the second term of Eq. (|B1[) 

E^^MIe-^Mff.l*). (B14) 
fc 

it is equal to a(t) for /3 p (0) = <5 p fc under Eq. (I3"5j) . Thus, 

^E^o^'-^'' 



(B15) 



i / dt' g' k e l( - u ' 
Jo 



Here, we will also multiply a factor exp(iCl't/2) to re- 
move the time dependent factor during the transforma- 
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tion. Therefore, 

4ReIp(t) = 4RcJ] A k (e, {0}| e - ,;ffrfft \g, l fc ) e l 



= 4Re[-i / dt'e ,!i ' i T4j[e-*'" i '] 

^0 j. 



2Q^ cos(Q' - u k )t - 1 

(Wfe + «) 2 fi'-Wfc 



2?rf / F(w,n / )G , 2(w)dw, (B16) 
Jo 



where 



= E 4 S + n? g ^ ( " _ Wfc) ' (B17) 

Here, we emphasize that the final result does not depend 
on the assumption that g k = 9 k - That's because it is the 
real parts of the second and third terms that contribute 
to the decay rate. 

For the fourth term in Eq. (|B1[) . we can also formally 
integrate a to have 



n = -i I dt'Y^g' k h{t')e 
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(B18) 



and substitute it into Eq. (136)) to yield 



dt'g'kY^dk'Pk'it' 



| e i(f2'-w fc /)t' e -i(0'-w fc )t 



rft ' e -«'(*-t')(^)2 e ic fc (t-t') ) 



where we have replaced 

P k ,(t')~P k ,(0) = 6 k , k . 
By one more iteration, we have 



(B19) 



(B20) 



/3 k (t) = 1 - / dt x / dt 2 (g' k ) 2 e^~ n '^-^. (B21) 
Jo Jo 



Thus, 



Y^A\ < 5 ,l fc | e -^«%,l fc > e * 



(B22) 



where 



is negligible since it is proportional to 
In total, 

xe 1 ^ 4 = Ci(l - 7 a ) - 2C 2 /,3 + ^^e- l(wfc ' 0,) * 



l-Z a -2/p + /a, 



(B23) 



where 



Z 5 =^^[e-^-^*-l], (B24) 



and we have dropped higher order terms of A\ . Then, 
the survival probability after one measurement is 

P(t) = \x(t)\ 2 

~ 1 - 2ReI a - 4ReI + 2ReI s 

~ exp(-2Re/ Q - 4Re/ /3 + 2Re/,5) (B25) 

with 
2Re/,5(i) =E' 



7k „ : „2 K - n')* 



sm 



/oo 
-OO 



, ./ v (w + n) 



2?rt / dwF(w,fi')G 3 (u;) 



sin — 



(B26) 



and 



G 3 M = E 



( w -tt') 2 a 



gi5{u - Wfc ). (B27) 



Straightforward, the survival probability after n con- 
tinuous measurements is 



2" _ „-m 



P(t) = \x(r)r = e 



(B28) 



where there are three contributions to the total decay 
rate 

R = Ei + R 2 + Ra, (B29) 

namely 

i?i = 2nReI a (r)/t 

/OO 
duF(w,il')Gi(u;), (B30) 
-oo 

R 2 = 4nRcIp(r)/t 

/OO 
dLoF(Lo i n')G 2 (uj), (B31) 
- OO 

R 3 = 2nReI s {T)/t 

/OO 
cZwF(w,fi')G 3 (w). (B32) 
-OO 

To conclude, the total decay rate is further simplified 



as 



/OO 
du>F(ui,Sl')G'(u)), 
-OO 



(B33) 
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where the measurement function is with the factor 

F(u,SV) = ^sinc 2 (^^r), (B34) 

ff ^1 . (3»-^ + 2^ fc )(O-r>0 
and the modified interacting spectrum Ty^k) — 1 t ^ ^ 2 • (cooj 

= Gi(w) + G 2 (w) + G 3 (w) 

= ^/( Wfe )^(a;- Wfe ) (B35) 
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